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We consider the nonequilibrium dynamics of the formation of a condensate in a spontaneously 
broken X(f>'^ scalar field theory, incorporating additive and multiplicative noise terms to study the 
role of fluctuation and dissipation. The corresponding stochastic Ginzburg-Landau-Langevin (GLL) 
equation is derived from the effective action, and solved on a (3 + l)-dimensional lattice. Particular 
attention is devoted to the renormalization of the stochastic GLL equation in order to obtain lattice- 
independent equilibrium results. 
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I. INTRODUCTION 

The importance of the investigation of phase transitions under extreme conditions was recognized long age 
efforts were devoted to develop a description within relativistic quantum field theory at finite temperature 
Among the variety of phenomena associated with the dynamics of phase transitions, phase ordering seems to be of 
particular importance for the understanding of the time scales governing the equilibration of systems driven out of 
equilibrium. The influence of the presence of a hot and dense medium on the dynamics of particles and fields is 
encoded "macroscopically" in attributes that enter stochastic evolution equations, usually in the form of dissipation 
and noise terms. Relevant time scales for different stages of phase conversion can depend dramatically on the details of 
these attributes. In particular, iiijiigh-energy heavy ion collisions, where one presumably forms a hot, dense, strongly 
interacting quark-gluon plasma ^5], chiral fields evolve under extreme conditions of temperature and energy density 
during the QCD phase transition. To have a clear understanding of data coming from BNL-RHIC, and especially of 
data that will be produced at CERN-LHC, one needs a realistic description of the hierarchy of scales associated with 
dissipation, noise and radiation, and also with the expansion and finite size of the system. 

Some of us have considered the effects of dissipation in the scenario of explosive spinodal decomposition (i.e., 
the rapid growth of unstable modes following a quench into the two-phase region of the phase diagram) for hadron 
production @, 0i H @] during the QCD chiral transition after a high-energy heavy ion collision in the simplest 
fashion llOl . Using a phenomenological Langevin description inspired by microscopic nonequilibrium field theory 
results T^. [Tsj. the time evolution of the order parameter in a chiral effective model [14] was investigated. Real- 
time (3 -I- l)-dimensional lattice simulations for the behavior of the inhomogeneous chiral fields were performed, and 
it was shown that the effects of dissipation could be dramatic in spite of the very conservative assumptions that were 
made. In fact, even if the system quickly reaches the unstable region there is still no guarantee that the unstable 
modes will grow fast enough. In particular, in the case of hadronization at RHIC time scales, the explosive scenario 
becomes rather constrained. More recently, analogous but even stronger effects were obtained in the case of the 
deconfining transition of SU{2) pure gauge theories using the same approach [Tsl. llGj. 

The framework for the dynamics was assumed to be given by the following standard GLL equation 

□ + 770 + V^/^(0)=C(x,i), (1.1) 

where is a real scalar field (the condensate {(f)) plays the role of the order parameter), Vg^j((/)) is the field derivative 
of a Ginzburg-Landau effective potential and rj, which can be seen as a response coefficient that defines time scales for 
the system and encodes the intensity of dissipation, is usually taken to be a function of temperature only, r/ = rj{T). 
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The function ^(x, t) represents a stochastic (noise) force, assumed Gaussian and white, so that (^(x, t)) = and 
(^(x, t)^(x', i')) = 2riT6{x — x')5{t — t'), in conformity with the fluctuation-dissipation theorem. 

However, a more complete field-theory description of nonequilibrium dissipative dynamics [ll[ shows that the 
complete form for the effective GLL equation of motion can lead to much more complicated scenarios than the one 
described by Eq. depending on the allowed interaction terms involving (/). In general, there will be nonlocal 

(non-Markovian) dissipation and colored noise, as well as the possibility of field-dependent (multiplicative) noise terms 

4> £, accompanied by density-dependent dissipation terms. In fact, on very general physical grounds, one expects 
that dissipation effects should depend on the local density (fy^cj) and, accordingly, the noise term should contain 
a multiplicative piece ~ cf). The typical term coming from fluctuations in the equation of motion for <j) will be a 
functional of the form 

T[(t){x)] ^ (l){x) j d^x'(j)^{x')Ki{x,x') + J d'^x'(j)ix')K2{x,x'), (1.2) 

where Ki{x, x') and K2{x, x') are nonlocal kernels expressed in terms of retarded Green's functions and whose explicit 
form depends on the detailed nature of the interactions involving (j). Explicit treatments for these nonlocal kernels 
show that under appropriate conditions one is justifled to express the effective equation of motion for in a local 
form 0, [H, [io, [2l|. The existence of these additional terms in the GLL equation will, of course, play an 
important role in the dynamics of the formation of condensates. For instance, it was recently shown that the effects of 
multiplicative noise are rather significant in the Kibblc-Zurek scenario of defect formation in one spatial dimension [2^ . 

In this paper we consider the nonequilibrium dynamics of the formation of a condensate in a spontaneously broken 
At/)** scalar field theory within an improved Langevin framework which includes the effects of multiplicative noise 
and density-dependent dissipation terms. The corresponding stochastic GLL equation is derived in detail from the 
effective action, generalizing the results of Ref. [TTj to the case of broken symmetry. The time evolution for the 
formation of the condensate, under the influence of additive and multiplicative noise terms, is solved numerically on 
a (3 -I- l)-dimensional lattice. Particular attention is devoted to the renormalization of the stochastic GLL equation 
in order to obtain lattice-independent equilibrium results. 

The paper is organized as follows. In Section II we review the perturbative approach of Ref. [Tl| . describe the 
effective theory that we consider, and present the computation of the effective action up to two loops in the Schwinger- 
Keldysh closed-time-path formalism, appropriate for a nonequilibrium situation. All results are extended to the case of 
spontaneously broken symmetry. In Sec. Ill we discuss the derivation of stochastic GLL equations with additive and 
multiplicative noises yielded by quantum corrections for a scalar A^* theory in the symmetric and broken symmetry 
phases. In Section IV the important issue regarding the proper lattice renormalization of the effective potential in 
order to achieve equilibrium solutions that are independent of lattice spacing is addressed. In Sec. V we show the 
results of our lattice simulations to study the behavior of the condensate. Section VI contains our conclusions and 
perspectives. One appendix is included in order to provide some technical details used in the paper. 



II. EFFECTIVE ACTION UP TO TWO LOOPS 



Here we compute the effective equation of motion obtained in Ref. [lli | for the symmetric phase, but instead of 
repeating the derivation in this reference, we make use of a more general method [l^, [13] that is valid for both 
symmetric and broken phases of the theory. The derivation can be easily extended to other field theoretic models. 

Let S[(t>\ be the classical action of a real scalar field 0, 



= j d^xC{4>), (2.1) 

with a Lagrangian density C{(j)) given by 

where TOq = rn? > in the symmetric phase and rn^ — —vn? < in the broken phase. In the presence of an external 
source J{x) that couples linearly to the field the classical action of the theory becomes 



S[(t>, J] = S[(t>] + j d^x J{x)(t>{x). (2.3) 
The corresponding generating functional Z [J] is given by 

Z[J] = exp '-W[J] ^ jT^4> exp S,[4>, J] , (2.4) 
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where the index c means that we perform our derivation in the context of the real-time Schwinger-Keldysh closed-time 
path formalism. In this formalism, the time integration is performed along a contour C with a branch C_|_ that goes 
from —oo to +00, and another branch, C_, that goes back to —00 [l^, [l^. The fields and sources on the branch C+ 
are denoted 0+ and J+, and those on the branch C_ are denoted 0_ and J_. Therefore, Sc\4>, J] in Eq. (|2.4p means 



Sc[4>, J] = S'[0+, J+\ - s[<j)-,j- 



(2.5) 



Our interest is the nonequilibrium effective action T[ip], which is a functional of a c-numher field Lp{x) containing 
all quantum effects. It is defined through the Legendre transformation 



where Lpix) is given by 



r[^] = w[j] - 

6W[J{x)] 



SJ{x) 



For convenience we have introduced the notation 

{J(j))^ = / d'^xJ{x)(f>{x). 



The external source can now be eliminated in favor of tp through 



J{x) 



S(p{x) 



(2.6) 



(2.7) 



(2.8) 



(2.9) 



The quantum effects encoded in ip can be made explicit by an expansion in powers of h, the loop expansion, which 
can be implemented in the following way [l^ . One starts by making the shift in the field variable (j) 



in the functional integral of Eq. p.4p , so that 



exp -r[. 



{x) (j){x) + (fix) 



V<p exp — 
h 



where we have used Eq. (|2.9p to ehminate J. Next, one expands 5'[(/> + Lp\ around = obtaining 

oo ^ 
n—l 

where S[ip] is the classical action and 



with 



Then we substitute this into Eq. (|2.1ip : 



(a;i) • • ■ d(j){xn) 



(2.10) 



(2.11) 



(2.12) 



(2.13) 



(2.14) 



A=0 



exp-r[(^] ^ I V(j) exp- 



2 {S2[^W)c 



2^(^nM</'")c-(flM(/.). 



n— 3 



(2.15) 



where we have made explicit the term quadratic (in 0) in the exponent of the functional integral, and have defined 

f [^] = r[^] - S[^] , (2.16) 
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and equivalently for the combination Ti[(p] — Si[(p]. 

Now, to make explicit the H dependence of the different terms in the expansion in Eq. (|2.15p we rescale (j) as 
(j) — > h^^'^cj) and expand r[(/3] and fi[(^] in powers of h, 



(2.17) 



n=l 



n=l 



to obtain the final result: 



expi^;i"-ir(")[v3] = V(j) expi 



1 °° f,n/2-l 



r!=3 



n=l 



(2.18) 



The functional T [ip] carries all the quantum corrections to the classical action S [tp] , without any approximation. In 
the following we consider the one- and two-loop contributions to the effective action. These are obtained by expanding 
both sides of Eq. (|2.18p in powers of h and equating terms of the same order. Specifically, the one-loop contribution 
is given by 



fW[(^] = -iln / exp ^ (52<^2)^ 



(2.19) 



Since the exponential is quadratic in 0, the functional integral in principle can be calculated exactly. However, in 
order to obtain an useful expression for deriving an effective GLL equation, we perform a perturbative calculation. 
This will be done in the next subsection. 

The two-loop contribution to the effective action is obtained keeping terms up to first order in h in the expansion 
of both sides of Eq. (|2.18p . The contributing terms are given by 



exp I 



(2.20) 



Expanding the ft-dependent exponentials on both sides of this equation, keeping 0{h) terms, and using Eq. (|2.19p for 
r*^^' [ip] into this last equation, one obtains for f [tp] the following expression: 



where the angular brackets {F{(j)) )c mean 



IV<PF{<P) exp(t52M0^), 
/^P0exp(f52M 0'), 



(2.21) 
(2.22) 



A. Perturbative expansion of r(i)[vp] 



So far, our treatment has been completely general in the sense that no differentiation between symmetric and 
broken phases has been done. However, for the implementation of a perturbative expansion it is necessary to make 
sure that the perturbation is performed around a stable field configuration. In order to ensure this situation and still 
keep treating both phases on the same footing we consider 'p{x) as 



ip{x) V + ipix) 



(2.23) 



where ^' is a x- independent field configuration. At tree level, in the broken phase where TOq = — < in Eq. ()2.2p . 
one has v = with = 6m^/A, whereas v — in the symmetric phase. In the following, though, we will not use 
the tree value of v. Instead, we determine it through the long-time behavior of the solution of the corresponding 
effective GLL equation. As we will show later, in the the symmetric phase v — Q, and in the broken phase v c^i ut^ for 
couplings not very large - which should be the case for the validity of the perturbative expansion. 
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For the evaluation of the one-loop expression for the effective action, Eq. (|2.19p . one needs the quadratic action 
{S2[ip]4>'^), where 5*2 M is defined in Eq. ^J^. We organize the perturbative expansion in a way that 

S2[^]=tS°{xi-X2)-XV"iip{xi))6^^\xi-x2), (2.24) 

where 

S°{xi - X2) = -i - m') - X2) , (2.25) 

and 

V'\v)^\v^^vV, (2.26) 



with 



= -m^ j^l^Xv^ , (2.27) 



Therefore, one can write for r(^^[(^] the expression 



r«M = 5M-iln / exp i (5V)^ 



X exp 



-i^A / d^xV"(ip{x))^^(x) 



(2.28) 



and the perturbative series is obtained from the expansion of the second exponential in powers of A. Specifically, 
defining 

Jvcl>exp^{S"cl>^)c, (2.29) 

and 

{F{<i,))l^^ jv<j, Ficj,) exp ^iSU\, (2.30) 

one obtains up to 0{X'^) 

r(i) [^] = SM-i\nZ~^jd^xV'\^{xm\x))l+i^jd^xV'\v{x))jdSV''{^{v)) [(</.(x)0(2/))O] ' . (2.31) 

To proceed one needs to make explicit the correlation functions ((/)(x)(/)(?/))2 on C+ and C-. There are four of 
such correlation functions, G'^^{x — x') with a,/3 = ±. Gq^ is the correlation hmction of (^{x) and (j){y) when both 
fields are on the branch C+ of the Schwinger-Keldysh contour, and therefore is the usual timc-ordcred propagator. 
Gq ^ is for both fields on C_, and gives rise to an anti-time-ordered correlation function. The G^^ and Gq refer to 
correlations of one field on one branch and one on the other branch. G^"*" is the usual physical (causal) propagator; the 
other three come as a consequence of the time contour, and are auxiliary quantities [29|. In the above, the subindex 
means that these correlation functions are calculated with the action 5°, defined in Eq. (|2.25p . All four G^^ can 
be expressed in terms of the usual advanced Gq (x, t) and retarded G'^ (x, t) correlation functions as 

G++(x,t) = Go>(x,t)0(i)-f G<(x,t)0(-t), 
Gq— (x,<) = Go>(x,t)0(-t) + G<(x,t)0(i), 
G+-(x,t) = Go<(x,t), 

Go-+(x,t) = Go>(x,t). (2.32) 
From now on, it is convenient to consider the spatial Fourier transforms of Gq'^(x, t), defined as 

Gf{^,t) = ^ I ^e^^-Gf{k,t). (2.33) 

The i on the r.h.s. of this equation is introduced for convenience, so that Gq^ (k, t) becomes real. In terms of these 
Fourier transforms, and making use of Eq. (|2.32p to express the G^^ in terms of Gq and Gq , the expression for r(i' [if] 
can be written as 
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^2 l-oo 

-i — / d^x d*x' 



d^q 



-ik- (x— x') 



(2^)^y (27r) 



k,0) 



V"i^+ix)) V"iip+ix')) G++(q, At) G++(q - k, Ai) 



V"{^+{x))V"{^^{x'))G+-{ci, At)G+-{q - k, At) - V"[^.{x)]V"[v+{x')] Gq +(q, A^Gq +(q - k, At) 



+ V"[v^{x)]V"{^^{x')) G„ -(q, At) Go -(q - k, At) 



(2.34) 



where At = t — t' . To proceed, one needs explicit expressions for G"^ . This wiU be discussed in Section Hill after we 
present the perturbative expansion for f [ip] in the foUowing. 



B. Perturbative expansion of r^^'[(^] 

Since {S2[^](f>^) in Eq. ()2.22p is quadratic, one can apply Wick's theorem and, after straightforward manipulations, 
one obtains the following terms that contribute to the expansion of r(2)[(^] up to 0{X^) 

f^'H^]=-^JfxV""{p) [(02(x))]%f(^Lg-su„M, (2.35) 

where 

f^cLg-sunM = ^^jd^xV"'{^{x)) 

X j^d^yV"'{^{v)) mx)<t>{y)),f . (2.36) 

and 

V"\p) = ip^v, V""{ip) = l. (2.37) 

To complete the 0{}?) expansion, one needs {(j){x)(j){y)) up to first order in A. This is given by 

(0(x)0(y)) = (0(x)0(2/)}o -txj^ d^z V"{^{z)) (0(x)0(z))o (0(2/)0(z))o • (2.38) 
Substituting this into the expression for r(2)[(^] in Eq. (jOS)) . one obtains 

f M --^J/'^ ^""(^) [('^'(^))o] ' + f (^'o„M + f lii„M , (2.39) 

where 

niL-.n.n[^]=^^Jf^V""ip) [{<jP{x)),] jd'x'V"{p{x')) mx)<j^{x')),f . (2.40) 

Using now the Fourier transforms of G"^, and following the same procedure as in the previous subsection of 
regrouping terms after using the expressions in Eq. (|2.32p . one obtains after a lengthy but straightforward calculation 

the following expression for rsnow-man [<y3] : 

fi^L-„.a„M = x/^'"*7 (^/ (i? ^""(^)Go++(k,0){F"[^+(.)] [G++(q,At)]^ 

^V"[ip.{x)] [G,+{ci,At)Y-V"[^+{x)] [G+-(q,Ai)] Vf"[^_(x)] [Gq "(q. At)] , (2.41) 
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where we used the fact that V""{lp+) = V""{(p-) = V""((p), in view of Eq. (|2.37p . In the same way, one obtains for 

f(2) r 1. 

setting-sun ir\ • 

X { V" [^+ {x)]V"' [^+ {x')] G++ (qi , At) G++ (q2 , At) G++ (qa , At) 

- V'"[ip+{x)]V'"[ip^ix')] G+- (qi, At) G+-(q2, At) G+-(q3, At) 

- V" {x)]V"' [^+ (x')] Go-+ (qi , At) Go + (q2 , At) Go + (qg , At) 

+ V'"[^-{x)]r"[^-{x')] Go -(qi, At) Gp— (q2, At) Gp— (qa, At)} . (2.42) 




Figure 1: Contributions to the finite-temperature effective action in the symmetric phase up to two loops and 0{X^). 




Figure 2: Contributions to the finite-temperature effective action, in the broken phase, up to two loops and 0{X^). 

In Figs. [Hand [2] we present the Feynman graphs representing the contributions up to two-loops and O(A^) to the 
finite-temperature effective action in the symmetric and broken phases, respectively. In these figures, each external 
leg means a factor of (p. It is easy to check that the graphs for the broken phase in Fig. [2] are those in Fig. [1] with 
each external leg shifted hy v - oi course the G"^'s are different in each set of diagrams, see Eq. (j2.25p . Notice also 
that if one were to obtain the graphs in the broken phase from those in the symmetric phase by shifting each external 
leg by I/, terms independent of the field variable ip would appear. These do not appear in the explicit and systematic 
calculation done above, but can be neglected because they do not contribute to the GLL equations. 



III. STOCHASTIC GLL EQUATIONS 



Here we derive an effective GLL equation. It has been shown in several previous works [ll|, |l7|, |3l|, |3^, l33| that 
dissipation is due to non-perturbative effects and cannot be derived from perturbation theory solely due to the 
appearance of secular terms 3J| which lead to the breakdown of perturbation theory well before dissipative effects 
can be observed. The issues involved in the derivation of the dissipation terms have been extensively discussed 
and examined, e.g in Refs. [l?] . [20l |. to where we refer the interested reader for further details. In particular, in 
these references it was shown that dissipation naturally emerges when the effects of quasi-particles in a thermal 
environment are taken into account, for instance, via the self-energy ^(q), its real part contributing to the dispersion 
relation E — i?(q) and its imaginary part to the width F = r(q) of the excitations. In this case, the correlation 
functions G^ and G^ can be written in momentum space as 



1 r 



G>{q,t) = -\ il + niEr))e 
G<(q,t) =G>(q,-t), 



-iErt 



n{EY)e 



iE*t 



(3.1) 



where Er = E — iF, and n{E) is the Bose-Einstein distribution. In the limit of F = 0, one recovers the usual 
expressions for a noninteracting Bose gas. The explicit expressions for the width F and the energy E can be found 
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in Refs. [Tl|, |T3|- From now on, we shall remove the subindex from C^^ to indicate that we are using dressed 
correlation functions. The fields ip^ and (y5_ can be expressed in terms of new variables and ip/^ defined as 



1 



■ifA + ^c, p>- ^ Pc 



(3.2) 



Using these into the expressions for f and r*^^\ one can express the result in terms of G++ only, and the effective 
action can be split into real and imaginary parts, F = Re F + i Im F. The explicit expressions for Re F and Im F for the 
symmetric phase are given in Ref. ll| . The corresponding expressions for the the broken phase can be obtained from 
those of the symmetric phase by the replacement of ipc{x) by = (fc{x) + v, and with the G^'^ that is appropriate 
for the broken phase. In view of this, instead of ipc we will use $ in the following, since when one takes v — Q, 
<i>(x) — > ^Pc{x) and the G"*""*" of the broken phase collapses to the symmetric phase expression when the concomitant 
change —m? is made. 

Given all this, the imaginary part of F[(/5Aj *&] can be associated to functional integrations over Gaussian fluctuation 
fields ^1 and ^2 by making use of a Hubbard-Stratonovich transformation [O, [l^l, leading to 



F[^A,<i>] = -t In J V^iP[^,] yi?6^[6] exp{z5e//bA,<i>,a,C2]}, 
with ScsiifA, 6] given by 

Scs[(pA,^,Ci,^2] = ReF[(^A,$]+ J d'^x(pAix)<i>ix)^i{x) + J d'^x ipA{x)^2ix) 
and P [^1] and P [^2] are Gaussian distribution functions given by 

PKi] - 7Vr'exp|-i J d^xd'x'^, (x) (^yRe [G++]l^}j 6 (x')^ , 

P[6] = iV^-iexpj-l J d^xd^x'i^ix) {^Re[G++]lJ^ , 

where N^^ and N^^ are normalization factors, and 



[G+ 



+12 ^ 



■ik- [ X— ) 



d^k 

(27r)3 " 

d^k (^"^i^-^') 



d\ 

3 . ,3 



G++(q,f-t')G++(q-k,i-0. 



.7=1 ^ ' 



(5(k - qi - q2 - qs) 



(3.3) 
(3.4) 

(3.5) 
(3.6) 

(3.7) 
(3.8) 



Given the distribution functions P [^1] and P one has that ^1 and ^2 have zero average and are correlated as 



(Ci(^)ei(^')) = YRe[G++]2^,,, 
(6(^)6(x')) = yRe[G++]^^,,. 



(3.9) 
(3.10) 



The fields ^1 and ^2 in (13. 4p act as fluctuation sources for the scalar field configuration if. While the field ^1 couples 
with both the response field fA and with the physical field $, leading to a multiplicative noise term in the equation of 
motion for the field ^2 gives origin to an additive noise term. As shown recently [2ll.[35j. both two-point correlation 
functions for ^1 and ^2 can be related to corresponding dissipation kernels appearing in the effective equation of 
motion for the scalar background field that can be derived from the effective action p.4p via (see e.g. [TTI. Isoj) 



'^■5'cff[yA,^,Cl,6] 



= 0. 



(3.11) 



This leads to a stochastic integro-differential equation given by 
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[□ + ml] ^[x) + A$3(^) ^ ^(^) ^ ^ $(^) ;C2[$] = + , (3.12) 

where /Cq ['i'] j ['^'l i a-nd /C2 are nonlocal memory kernels given by 

^ , ,3.3, 

/Ci[<i>] = y 1 d'x' j' dt'$(x',t')Im[G++]^ ^, , (3.14) 
^2[1>] = yl d^-T'l* dt'$2(x',i')Im[G++]^^^, . (3.15) 

We reiterate that this equation is valid for both symmetric and broken phases of the theory. 

One problem one has to face here is the renormalization of ultraviolet divergences. Although the renormalization 
of ^efi is defined by the usual introduction of counterterms in the initial Lagrangian, the nonlocal memory kernels 
involve additional complications. In particular, the numerical simulation of nonlocal equations like Eq. (j3.12p is a 
very difficult problem. However, the situation becomes simpler after a series of physically motivated approximations 
related to spatial and temporal non-localities of the memory kernels. Considering only contributions which are slowly- 
varying in space and time (this is a valid assumption for systems near equilibrium, when $ is not expected to change 
considerably with time - quasi- adiabatic approximation [3ll[3^1. and keeping contributions up to order A^, Eq. p.l2p 
can be put into a GLL form ^llj 

(□-Hm2,)$(x,t) + y$3(^^^)^^^$2(^^^)^(^^^) ^ ^3^g) 

where 771 is a temperature-dependent dissipation coefficient, At and ttit (see below) are the renormalized finite- 
temperature coupling constant and mass, respectively. As explained in [ll|, the additive noise term, ^2, appearing 
in Eq. p.l2p drops out at C'(A^), as well as the corresponding dissipation term, ^ rj2- Despite the approximations 
leading to Eq. (j3.16p may imply in a number of issues that need to be accessed for the validity of these approximations 
(see e.g. jlTj). in our numerical simulations we adopted Eq. (I3.16P so to make the analysis simpler and so assume the 
dissipation coefficient rji as giving an order of estimate for the magnitude of dissipation in the system. An analysis of 
the validity of Eq. (|3.16p based on the basic assumptions leading to it will be checked in Sec. V.B. 

The parameters in Eq. (|3.16p are given as follows. For the symmetric phase (ttzq = > 0) one has in the 
high-temperature limit T/m ^ 1 and within the validity of perturbation theory in the thermal theory, XT/m <C 1, 
that niT and At are given by 

ml ~ m^-fA— , 3.17) 
24 



At — A 



..(3^)+0(a'^), (3.18, 



167r2 

and the dissipation coefficient 771, also in the high-temperature limit, given by [Til. [itI [3l| 

T»mT 96 / T A , , 

ryi c."— In 2— . 3.19 
TTJ \ mx J 

For the broken phase, ttiq = — < and the expression for m^ changes in view of the possibility of symmetry 
restoration at sufficiently high temperatures. Specifically, the dominant terms are now given by 

9 T»mT 9 AT^ 

mi ~ -m^ + — - -(-•■•. 3.20) 

24 

2 — 0/1^t2 , 



From this, one has a critical temperature given by — 24m /A 
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IV. LATTICE RENORMALIZATION 

In our following numerical study, we consider an extended GLL equation, incorporating additive and multiplicative 
noise terms. The time evolution of the field ip(x, t) at each point in space, which will determine the approach of the 
condensate {lp) to its equilibrium value will be dictated by the following equation: 

□ (^+(r/i if" +r^2)^ + V,ff{v)^iiM <^ + 6(x,t) , (4.1) 

where -qi and 772, which can be seen as response coefficients that define time scales for the system and encode the 
intensity of dissipation, will be taken to be functions of temperature only, 77^ = VaiT) (a — 1,2). The functions ^i(x, t) 
and ^2(x, t) represent stochastic (noise) forces, assumed Gaussian and white. 

A. The need for lattice counterterms 

Analytic solutions of Eq. (j4.1[) are achievable only in very special situations, like in a linear approximation to 
Veff{(p), usually valid only at short times. Complete solutions describing the evolution of the system to equilibrium 
can be obtained only through extensive numerical simulations. In general, numerical simulations are performed on a 
discrete spatial lattice of finite length under periodic boundary conditions. However, in performing lattice simulations 
of Eq. (j4.ip . one should be careful in preserving the lattice-spacing independence of the results, especially when one 
is concerned with the behavior of the system in the continuum limit. The equilibrium probability distribution for 
the field configurations ip that are solutions of Eq. (|4.ip is Peq[y^] = e^^I'^', where S[ip\ is the Euclidean action. The 
corresponding partition function is given by the path integral 

Z[ip] = J'^f e~^''^' • (4.2) 

The calculation of expectation values and correlation functions of ip with this partition function leads to ultraviolet 
divergences. In the presence of thermal noise, short and long wavelength modes are mixed during the dynamics, 
yielding an unphysical lattice spacing sensitivity. The issue of obtaining robust results, as well as the correct ultraviolet 
behavior, in performing Langevin dynamics was discussed by several authors [36l. IstI. [ssj . The problem, which is not 
a priori evident in the Langevin formulation, is related to the well-known Rayleigh- Jeans ultraviolet catastrophe in 
classical field theory [3§|. The dynamics dictated by Eq. (|4.1I) is classical, and is ill-defined for large momenta. 
These a priori lattice divergences can be eliminated by renormalizing the potential Vg// through the addition of 
appropriate counterterms (notice that these divergences are completely unrelated to the usual ones of the quantum 
theory). Since the divergent terms are all perturbative, one can identify the appropriate diagrams, subtract their result 
computed within the classical theory, then add the equivalent terms computed within the quantum theory. In our case, 
this corresponds to substituting the Maxwell-Boltzmann weight by the Bose-Einstein distribution in the momentum 
integrals for the appropriate divergent terms. Since the theory in three-dimensions is super-renormalizable, only a 
mass renormalization is required - only the one-loop tadpole and two-loop setting-sun diagrams are divergent. We 
only require, then, a renormalization of the mass parameter tuq in Ve// ((/?). Using such a renormalized potential in 
Eq. (|4.ip leads to equilibrium solutions (p that are independent of the lattice spacing as we are going to explicitly show 
below. In practice, this lattice renormalization procedure corresponds to adding finite-temperature counterterms to 
the original potential Vef / {<p) , which guarantees the correct short- wavelength behavior of the discrete theory as was 
originally shown by the authors of Ref. ^36^ within the framework of dimensional reduction in a different context. 
Furthermore, it assures that the system will evolve to the correct quantum state at large times. 

In order to define the notation and the procedure to calculate the loop expansion of the effective potential [l^l , we 
first consider the theory in the continuum and switch to the lattice when calculating the loop integrals. 

B. The classical field effective potential in the continuum: loop corrections 

The calculation of the classical three-dimensional loop corrections to the potential starts as usual by introducing a 
constant field ip through the transformation 



(j)^(j) + ip, 



(4.3) 
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and considering the action 



S[ip; 0] = S[ip + (f>]- S[^] - £xd) 



dS[4>] 



del) 



(4.4) 



Next, the quadratic terms in S are collected in a "free" action 52, and the remaining terms in an interacting action, 
Si. Then, the effective classical field potential is defined by the expression 



(4.5) 



where V is the three-dimensional volume. From this expression, one obtains for V, 



effl 



^ln(e-^^[-^l). 



(4.6) 



where 



-SiVp,<t>]\ — 



(4.7) 



The final step is the determination of ip through 



dV, 



effi 



dtp 



(4.8) 



The action 5*2 [<^; 4>] for bare potential under consideration is given by 



— TO 

2 



2 J,2 



(4.9) 



where 



2 2 , \ 2 



(4.10) 



Since ip is constant, the first functional integral in Eq. (|4.6p can be easily performed in momentum space, with the 
result 



(3V 



In / V(j)e 



-S2[v;'t>] 



T 



lnG-i[(p;fc2], 



(4.11) 



where G is the inverse of the three-dimensional (classical field) propagator 



fc2 



(4.12) 



and = / . The calculation of the setting-sun diagram requires Si, 



Si[^-A]^(i I d^x (-1^03^1a0^) , 



(4.13) 
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with K = — iAiy9. Expansion of e ^' in Eq. (|4.6p gives 



1 - t2 / 1 \ 2 



where 



6V 



G[ip;k']G[ip;q']G[ip;ik + qy]. (4.15) 



' /c g 

Now, the divergent part of Vef / [<^] is obtained from 



tadpole+ 2 

setting— sun 



XT X^T^ 

Imv[^] —Hd^v[^], (4.16) 



I[ip] = / G[^;fc2], (4.17) 



with 



where Idiv[^] and i?dij,[iy9] represent the divergent parts of I[ip\ and H[ip]. Notice that the derivatives of H[ip] with 
respect to ip lead to finite integrals and hence are irrelevant here. The evaluation of the divergent parts of / and H 
requires a regularization scheme. Since we simulate our GLL equations on a cubic lattice, we shall evaluate these 
divergent parts by calculating the effective potential on the lattice. 

C. The counterterms for the classical field potential on the lattice 

Here we consider the theory on a cubic lattice of volume V — L'^, with L = Na, where a is the lattice spacing and 
iV is the number of lattice spacings. The coordinates Xi of the lattice sites are such that < < a{N — 1). The 
Laplacian on the lattice, V? = A, is defined as 



1 ^ 

A(t){x) = —^^4>{x + ai)-2(l){x)+(j){x~ai) , (4.18) 

^ 4=1 



where i is the unit cartesian vector indicating the three orthogonal directions of the square lattice. We also impose 
periodic boundary conditions (PBC) on the fields, 

(j){x + aNi) = (j){x) , (4.19) 
and define the Fourier transform /(fc) of a function f{x) on the lattice as 

/(fc) = a3^e-*V(x). (4.20) 
Because of the PBC, the allowed lattice momenta form the Brillouin zone B 

2tt 

h = —n,, n, = 0,1,2,- •■ ,7V- 1. (4.21) 
aN 
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The inverse transform is given by 



(4.22) 



fc eB 



and the momentum summations over the Brillouin zone will be denoted by 



V 



k eB 



Let us now consider the derivation of the effective potential on the lattice. The lattice action is given by 



(4.23) 



(4.24) 



where V[(p] is the same as before. The derivation of the effective potential follows the same procedure as in the 
continuum, leading to expressions for the setting-sun and tadpole graphs as in Eqs. (|4.15p and (|4.17p . but with J^, 
given by the sum over lattice momenta as indicated in Eq. (|4.23p . What remains to be determined is the lattice 
propagator corresponding to G. The lattice propagator Giatt can be obtained from the quadratic action S2[^',4'] on 
the lattice 



—m 
2 



(4.25) 



with vn? given by Eq. (|4.10p . The lattice propagator Giatt[VT^Ty] is the inverse of (—A + m^), i.e. 



(4.26) 



where 5x.y is the Kronecker delta. Since 1^9 is a constant field, the solution of this equation can be obtained using the 
Fourier transform of Giatt['~P', x, y\, defined as 



Giatt[v]x,y\^ / e'^'''GiattYp]k 



Substituting this into Eq. (|4.26p . one obtains 



4 d(ip;ni,n2,n3) ' 



where 



(4.27) 



(4.28) 



d{(p;ni,n2,n3) = — 



2a ^ El (1 ~ cos akj) - 
E sin^{Tm,/N) + {am/2f 



(4.29) 



One can then write I[ip\ for the tadpole diagram as 



1 ^ 



(4.30) 
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and the double sum in H[(p] for the setting-sun diagram as 



E 



^ d{(p;ni,n2,n3)d{(p;mi, 1712,1713) 

1 

m + TOi, 712 + "^2, '^a + 'tis) 



(4.31) 



The divergent parts of the sums above can be isolated in the limits of TV — > cx) and a — > 0. For example, the tadpole 
sum in the limit of — > 00 can be converted into an integral [36j . 



/M = 



1 



4a7r-^ 7q 
1 



1 



sin + (am/2) 



4a Jo 



e-"/2/o(a/2) 



where /q is the modified Bessel function. In the limit of a ^ 0, the divergent part of is given by [36 



(4.32) 



where S is a constant defined by 



+7r/2 



-7r/2 



Ei sir 



3.1759 



(4.33) 



(4.34) 



It is important to remark here that for T < Tc it is crucial the use of Eq. (j4.8p so that rr? < in the above expressions. 
In a more demanding numerical computation, Jidiv \^\ can also be isolated, with the result [36l | 



167r2 









+ c 





(4.35) 



where M is the renormalization scale and C ~ 0.09 is another constant appearing in the integrals being evaluated 
numerically [3^. 

The renormalized mass in the classical field perturbative loop expansion is obtained as in the case of quantum field 
theory, i.e. by subtracting the divergent parts of these graphs as indicated in Eq. (|4.16p : 



1 



2 i2 

m (p 



1 



2 i2 



where the mass counterterm is given by 



6m^ 



. f 1 M [ 1 

= -z- idiv m — ^ Hdiv m 



2 "■^'"^^^ 6 

Therefore, we add the following finite-temperature counterterms to our original potential: 



Svr a 96tt^ 



V aM / ^ 



(4.36) 



(4.37) 



(4.38) 



Note from the above equation that the dependence on the mass scale M of the lattice conterterm is only logarithmic. 
So it turns out that results are only very weak dependent on changes of the scale M . Of course, any change in 
the renormalization scale M can be compensated by corresponding changes in the renormalized parameters (which 
here requires only a change in the renormalized mass, as it is clear from Eq. (|4.38|) ). and as also expected from the 
renormalization group theory. 
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V. NUMERICAL RESULTS 



In this Section we present results of extensive numerical simulations for the GLL equation obtained in three spatial 
dimensions. The numerical methods used to do so are described in the Appendix. 

Nowadays the understanding of the relaxation of a field to its equilibrium configuration is of great importance, once 
there are many natural phenomena that occur out of thermodynamic equilibrium. In order to analyze this behavior, 
we start by studying the dependence of the solutions on the lattice spacing a = L/N and then, by introducing the 
lattice counterterms discussed and derived in the previous section, we show how lattice-independent results can be 
obtained in Langevin simulations using both standard additive noise and dissipation terms and also in its generalized 
form, which includes field-dependent (multiplicative) noise and density-dependent dissipation. This is particularly 
important, since as we have shown in Sec. Ill, the expected effective equations of motion for background scalar fields 
result to be in general of the generalized Langevin form. 

All results that will be presented here refer to the time dependence of the volume average of the order parameter, 
defined as 



(<^(a:,j/,z,t)) = ^^^r,fc, (5.1) 

ijk 

where t/)"^ is the average over a large number Nr of independent noise realizations. 



r=l 

In all our numerical Langevin simulations, we consider Nj. between 20 and 100. As usual, lattices with larger values 
of L require relatively less realizations over the noise. We have considered and tested different lattice sizes to ensure 
the robustness of all numerical results. 



A. The problem of lattice dependence in the generalized GLL approach 

As discussed in the previous section, the simulation of equations with noise, being classical by nature, leads to the 
appearance of Rayleigh-Jeans ultraviolet divergences at long times when simulating the equation on a discrete lattice. 
These divergences manifest themselves in the form of lattice-spacing dependence of the equilibrium solutions. We can 
show that by just considering the easiest nonequilibrium evolution, which is the one of relaxation to the equilibrium 
vacuum state with initial conditions away from it. In Figs. [2]and|l]we show the corresponding dynamics for the scalar 
field expectation value, {4>{x,y,z,t)) defined in Eqs. (|5.ip and (|5.2p . in the broken symmetry case, T < Tc, with the 
critical temperature Tc extracted from the finite-temperature effective potential. 

The initial state for the field in the simulations for the broken phase was taken around the inflexion (or spinodal) 
point of the finite-temperature (Ginzburg-Landau) potential, 0i„fl, defined by 



d(lP- 



= 0. (5.3) 



In Fig. [3] we show results for the standard Langevin equation (jl.ip . while Fig. H] displays results for the generalized 
case. The parameter values considered for the temperature and dissipation terms 771 and 772, in units of the scalar field 
mass, were T/m = mrji = 1^2 I'm — 1, while the quartic coupling constant was chosen to be A = 0.25. The scale M 
is taken as M/m = 1. These values suffice for our purposes of just demonstrating the lattice dependence problem in 
Langevin simulations. In Figs. [3]and[l]the number of lattice points and the time stepsize were kept constant, = 64 
and 5t = 0.01 (in units of the scalar field mass), respectively, while the lattice spacing, mSx = a — L/N, was varied. 

It is clear that the solutions shown in Figs. [3] and |4] are not stable as the lattice spacing is modified. As discussed 
previously, this problem can be traced to the fact that the equilibrium value of the quantity ipijk gives the classical 
average 



16 



6 



5 



1-4 



3 



^0 5 10 15 

m t 

Figure 3: Solution of the standard GLL equation using the leap frog algorithm for different lattice spacings. 
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Figure 4: Solution of the GLL equation with both additive and multiplicative noise and dissipation terms using the leap frog 
algorithm for different lattice spacings. The other parameters are the same as in Fig. O 

an ultraviolet divergent quantity. Thus, stable equilibrium solutions of the GLL equation, i.e., solutions not sensitive 
to lattice spacing, can only be obtained by the introduction of the appropriate counterterms in the effective potential 
in order to eliminate these divergences. These counterterms are the ones derived in the previous section and given by 
Eq. gJll). 

In Figs. [5] and [6] we present the results of the simulations including the counterterms. As we can see, equilibrium 
solutions that are independent of lattice spacing are obtained. 

For the standard GLL equation this was shown extensively in a series of previous papers [STj , but we are not aware 
of the same demonstration for the case including multiplicative noise terms. From the results shown in Figs. [5] and [6] 
we can also immediately draw a couple of interesting conclusions. The first is that even though the counterterms used 
were calculated with an equilibrium partition function, thus ensuring lattice-spacing independence only in equilibrium 
(long times) situations, the results for short times show only small lattice-spacing dependence. Second, we can notice 
that the relaxation time scales to the vacuum state is about the same in all cases. Another observation we can make 
based on the results shown in Figs. [5] and [6] is that the dynamics with multiplicative noise and dissipation terms can 
be quite different from that driven by additive noise only. In particular, notice from Fig. [Hlthat the relaxation time 
to the vacuum state with the generalized GLL equation is much longer than the one with just additive noise, even 
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Figure 5: Solution of the standard GLL equation using the leap frog algorithm for different lattice spacings and including the 
renormalization counterterms. 
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Figure 6: Solution of the GLL equation with both additive and multiplicative noise, and dissipation terms and including the 
renormalization counterterms. 

though the magnitude of the dissipation coefficients (in the dimensionless units in terms of the field mass m) are of 
order one. Also, no oscillation of {4>) is seen in Fig. [6l When we consider the typical magnitude of the multiplicative 
dissipation coefficient obtained within the local approximation, Eq. (I3.19p . which is about two orders of magnitude 
larger than the one used in the example of Fig. [6l we expect a much longer relaxation time for the scalar condensate. 
This is typical of a regime of strong overdamping, implicit in the derivation of the local equation (j3.16p p7| . As we 
are going to see next, the study of Eq. (j3.16p on the lattice will provide us a unique way of accessing the limit of 
validity of the approximations used in order to arrive at that equation, thus generalizing to a much more realistic 
situation the simplified studies performed in Ref. p^j . 



Simulation results and tests of validity 



Finally, we present results of simulations for the GLL in the symmetric and broken cases. Let us first study 
the broken case to complement the results presented for the phenomenological parameters considered e.g. in Fig. [6l 
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Considering Eq. (|3.16p . with dissipation coefficient given by Eq. (|3.19p . we show in Fig. [7] results obtained for different 
values of temperature T < for a cubic lattice with = 64 and L — 32. The stepsize in time is 5t — 0.01 (as before, 
we normalize all quantities by the scalar field mass m) and in all results to be presented henceforth, the proper lattice 
counterterms to the effective potential, Eq. (|4.38p . were taken into account. The initial value for the field in all cases 
was around the inflexion (spinodal) point, 4>iri&, Eq. (|5.3[) . at each respective temperature. 
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Figure 7: Solutions for different values of temperature of the microscopically derived effective equation of motion in the broken 
phase. 

The curves displayed in Fig. [7] show the large relaxation time to achieve the equilibrium vacuum state due to the 
large magnitude of the multiplicative dissipation term 771, as anticipated in the previous subsection. We can also see 
from the results of Fig. [7] that the dynamics for the scalar field condensate is mostly overdamped all the way to the 
equilibrium state, except by some very small transient nonequilibrium evolution around the initial time. This seems 
to support the adiabatic approximation (of slowly varying field) used to derive the local effective equation of motion 
for the condensate field and used to obtain e.g. Eq. (|3.16p from the full nonlocal result, Eq. p.l2p . It shows that 
the lattice evolution of the nonhomogeneous field (/>(x, t) can be much less restrictive with respect to the adiabatic 
approximation than the one for the homogeneous field approximation considered in Ref . T?] . In this reference it was 
shown that the dynamics for a homogeneous field could only be overdamped in the high-temperature, high-amplitude 
regime, and when the scalar field is coupled to a large number of other scalar fields. The results shown here seem to 
indicate that these restrictions are not all necessary for an overdamped dynamical regime. As for the other conditions 
used to derive Eq. (|3.12p . and then the local approximation (|3.16p . are concerned, e.g. the amplitude expansion for 
the field, and the high-temperature approximation used in the derivation of the dissipation coefficnt (j3.19p . these are 
clearly valid for the parameters considered in our case to obtain Fig. [T) Other values of coupling constant A and 
initial field value can easily be found to satisfy these conditions as well. 

As for the basic approximation used in all calculations shown in Sec. II, i.e., the loop expansion, our lattice 
simulations also offer a simple test to check its validity. Recall that the loop expansion is valid provided that the 
(vacuum and thermal) fluctuations at some given temperature are sufficiently small, in which case the classical field 
configuration used to derive the functional partition function is close enough to the solution that extremizes the 
classical action. This is not so if the fluctuations have large amplitudes. In Ref. [iH it was shown, for the broken 
potential, that these fluctuations should have an amplitude not larger than the inflexion point of effective potential and 
in that reference, fluctuations were modeled as subcritical bubbles (see also [i^l) in order to compute their amplitude, 
or actually their root mean square (RMS) amplitude. The equivalent quantity that can be considered here is the 
RMS variance of the field on the lattice, defined as 



-/-RMS = VW)^W , 

where the averages are on the lattice and taken analogously as in Eqs. (|5.ip and (|5.2p . 



(5.5) 



19 




400 600 

m t 



1000 



Figure 8: The RMS amplitude of fluctuations (normalized by the equilibrium (^inti) as a function of time. 



In Fig. [5] we shown 0RMs(i) for the four temperatures considered in Fig. [71 We see that 0rms > ^infl already for 
T > O.TTc, beyond which thermal fluctuations become too large and the loop expansion used to derive Eq. (|3.12p is 
no longer valid. 

It is important to notice that the thermalization of the field in our simulations is a rather fast process, happening 
typically in a very short time, ithorm ~ 0.5/m, which is mostly due to the large dissipation coefficient p.l9p coming 
from our calculations and subsequent approximations. During the transient initial time given by tthorm, where it is 
mostly out of equilibrium, the background field quickly reaches a thermal equilibrium stationary state. Thus, the 
thermalization time in a Langevin simulation is a rather quick process, even though the system is still evolving in 
time. The actual equilibration time for the field in its vacuum state happens in a much longer time, as observed in 
our simulations. 

Next, we consider the symmetric phase (with positive mass term in the potential). Results for the symmetric 
case are presented in Fig. ^ For this case, we considered parameter values satisfying the perturbative expansion 
(m^ 3> A(/)^/2, XT/m < 1) and high-temperature approximation, given by A = 0.01, T/m = 5, and initial field 



amplitudes 4>o — 5m and (/)q — 10m, such that we could compare to similar values used in Ref. [17[ in order to probe 
the overdamped dynamical regime for the scalar field. The lattice parameters are the same, N = 64, L = 32, as well 
as the stepsize in time, St = 0.01. 

Notice from Fig. [9] that the dynamics is again mostly overdamped, with the relaxation time for the field increasing 
with the initial field amplitude, thus seeming to support the results of Ref. [13] that showed that the overdamped 
regime is more favorable the larger the field amplitude is. As for the broken case, we can again probe the validity 
of the loop expansion used to derive the effective evolution equation for the field. Since in this case the potential 
has always a positive definite curvature, we use the conservative criterion that the perturbative expansion is valid for 
fluctuations dominated by the quadratic term of the effective potential, i.e.. 



Ofluct 



< 



(5.6) 



where 



= W12- 



(5.7) 



with m^ and At given by Eqs. (|3.17p and (|3.18[) . respectively. In Fig. [10] we show the RMS amplitude of fluctuations 
of the fleld on the lattice for the two cases shown in Fig. [5) Since, from Eq. (|5.7p and the parameters used, (/)max ~ 34m, 
both cases shown in Fig. [Sj comfortably satisfy the criterion (|5.6p for the validity of the loop expansion. 
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Figure 9; Solutions of the microscopically derived effective equation of motion in the symmetric phase for two different values 
of initial field amplitude. 
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Figure 10: The RMS amplitude of fluctuations of the field on the lattice for the two values of initial field amplitudes shown in 
Fig. [9] for the symmetric case. 

VI. CONCLUSIONS 

In this work we have studied several important aspects regarding the dynamics of a scalar field background con- 
figuration in both symmetric and broken symmetry phases. It has been long recognized that the effective evolution 
equation for the field can be of a complicated form. From the analogy with standard Langevin equations for the study 
of the approach to equilibrium, the microscopically derived effective evolution equation allows for the presence of 
similar additive noise and dissipation terms, but also for multiplicative (field-dependent) noise and dissipation terms. 
Although equations of motion of the standard Langevin form (with only additive noise and dissipation) have been ex- 
tensively studied in the literature, its generalized form, which includes the multiplicative noise and dissipation terms, 
still demands extensive studies. Here we have performed a number of numerical simulations with these equations 
on a cubic lattice. We have also called the attention to another issue frequently overlooked in the literature: the 
necessity of adding lattice renormalization countertcrms to cancel Rayleigh- Jeans divergences of the corresponding 
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classical theory, in order to produce sensible results from Langevin simulations. We have shown that the same lattice 
counterterms that are required for the standard Langevin simulations also work for the generalized Langevin equa- 
tions, producing lattice-independent equilibrium quantities, and also minimizing the dependence of the dynamics on 
the lattice parameters. 

In our studies of the generalized Langevin equation for both the broken and symmetric cases, with a typical value 
of dissipation coefficient as produced from a local and high-temperature approximation for the effective equation of 
motion for the background scalar field, we have shown that the dynamics produced is overdamped. Moreover, we were 
able to test some of the approximations used to derive these equations, such as the validity of the loop expansion. 
This allowed us to explicitly test some of the hypotheses and results obtained in Ref. , studied there only for the 
homogeneous field approximation. 

One interesting and important problem that still remains, though beyond the objectives set for the present work, 
is the study of the validity of the approximation of transforming the complicated nonlocal (non-markovian) equations 
of motion obtained through a microscopic derivation (via the Schwinger-Keldysh real-time formalism for the effective 
action) into the local form used in the simulations performed for this study. This is a notoriously difficult problem 
due to the oscillatory nature of the nonlocal kernels appearing in the full effective equation of motion, which leads to 
uncontrollable numerical behavior in simulations. We are currently studying viable numerical strategies to deal with 
this problem and to fully access the local approximation for the equations, as well as the differences in its dynamics 
as compared to the nonlocal equations. Results in this direction will be presented in a future publication. 
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Appendix A: THE GLL EQUATION ON THE LATTICE 



The equation to be solved on the lattice is 



_ V^^ + r;i</>2^ +^,^+ (0) = + 6 . (Al) 

We solve the GLL equation on the lattice with periodic boundary conditions (PBC). When we insert the system in 
a box, (j) acquires a discrete form 4>2jk^ where t = nAt with n = 0, 1, 2, . . ., x = ia, y = ja and z = ka, being a the 
lattice spacing a= Using PBC we have 



VN+ljk 
^N+lN+lk 
^N+lN+lN+1 
€,0 



Firstly we write the Laplacian 
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f 1 ' 
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■'iOO = 
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i^iNk ' 
4^NjN 



(A2) 



-Ijk ~ '^ijk 
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(A3) 
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For simplicity, we introduce the compact notation 



ijk 



(A4) 



In addition, we divide the time in n steps, t = nAt, with n = 0, 1, 2, . . .. With respect to the discretization of the 
time derivatives, we use the leapfrog approximation method, where the algorithm is defined by the following iteration 
scheme: 



dt 



^71+1/2 



At 



At 



9ji+l/2 ~ <Pji-l/2 



(A5) 



Now we rewrite ()Al[) in terms of discrete field (pf^f. 



Using the discretized quantities derived above we obtain 



^"ifcCl + $2 



(A6) 
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(A7) 
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The equation that we solve by iteration is 



m{€jk) +772 



;n+l/2 
hjk 



At 



At 



??2 



'--2 



mi^.kY + m A<U"-i/^ + (L</.)" At 



V^s (</>Sfe)At-fC,fc6At + 6At 



In a compact notation we have 



(A8) 
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;n+l/2 
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(AlO) 



where 
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m {<KokT + m 



At 
At 



(All) 



To update the field we make 
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(/.I;+l = c,, + A^<^:;+^/^ (A12) 

Finally, on the lattice the noise terms are modelled to satisfy the discretized fluctuation-dissipation relation {i — 1,2 
denote one of the two noise and respective dissipation terms appearing in the generalized GLL equation (j4.f p ) 

{^^,n^J,n') = 2^,T S , J n ,n' / {a^ At) , (A13) 

so that its amplitude that can be written as 



where Gi^n is obtained from a zero-mean unit-variance Gaussian. 

Finally, since we are dealing with multiplicative noise in a Langevin equation, there is a well-known ambiguity 
in the discretization of time [43| . This ambiguity is manifest whenever there are multiplicative noises and they are 
delta-correlated in time (i.e., they are Markovian, like in the approximation we have adopted in our simulations), 
and leads to the two most common types of discretization procedures, the one by Ito and the one by Stratonovich 
[4^ . The leapfrog scheme shown in Eq. (|A10|) corresponds to Ito's prescription. The ambiguity can be seen when 
one discretizes the (Markovian) multiplicative noise term. Notice from Eq. (IA14P that the noise terms in Eq. (jAlOp 
are actually of order (At)^/^, instead of order At as the remaining terms. So, the multiplicative noise term needs 
to be re-expanded up the to next order in At. This can be performed when one writes the discretized form for the 
multiplicative noise term using the Riemann formula 

/t+At 
(j){x,t')£,i{x,t')dt' = [{I - a)(f>{x,t) + a(l){yi,t + At)] xii.x,t) , (A15) 

where < a < 1, with the Stratonovich prescription corresponding to a = 1/2, while Ito's corresponds to for a = 
p3 |. In (|A15p xi(x, i) is a new Gaussian stochastic process described by (in discretized form) 

iXl.nXl.n') = 2inT5n^n'At/a'^ . (Al6) 

Using Eq. (jAlSp back in the leapfrog equation (jAlOp and re-expanding it to the next order in the multiplicative 
noise term one finds that in the Stratonovich interpretation there is an additional correction term to the leapfrog 
equation (jAlOp of order a(j)'^~^}^'^ xiAtQ /'E? , which is already an order (At)^/^ higher than the remaining terms in 
that equation. We have explicitly checked in all our simulations that this is a negligible correction for our results, and 
so we have adopted the Ito's prescription, i.e. Eq. (jAlOp . Notice that this might not be the case would we be working 
with first order in time derivative equations, in which case large corrections due to the difference between Ito's and 



Stratonovich's interpretations can arise [45l|1. For a recent discussion on different prescriptions in the context of the 



relativistic Brownian motion, see e.g. Ref. (46|. 
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